Electronic continuum states and far infrared absorption of InAs/GaAs quantum dots 
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The electronic continuum states of InAs/GaAs semiconductor quantum dots embedded in a 
GaAs/AlAs superlattice are theoretically investigated and the far/mid infrared absorption spec- 
tra are calculated for a variety of structures and polarizations. The effect of a strong magnetic 
field applied parallel to the growth direction is also investigated. We predict that the flatness of 
the InAs/GaAs dots leads to a mid-infrared absorption which is almost insensitive to the magnetic 
field, in spite of the reorganization of the continuum into series of quasi-Landau states. We also 
predict that it is possible to design InAs/GaAs photoconductors which display very strong in-planc 
absorption. 



PACS numbers: 73.21. La, 73.21.Cd, 78.67.Hc 
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INTRODUCTION 



It is by now well established that the lowest lying eigen- 
states of the InAs/GaAs quantum dots (QDs) are bound 
states if the dots are not too small. At higher energies 
one finds continuum states which arc cither bulk-like (i.e. 
display extended wave functions in the three directions of 
space) or are bound along the growth axis, the wetting 
layer statesi. The continuum states are of paramount 
importance to the control of useful electronic properties 
associated with the bound states (e.g. lasing action or 
far infrared photoresponse). This is because their density 
of states are considerably larger than that of the bound 
states of the dots since the QDs are most of the time 
very diluted. Despite their importance, there have been 
relatively few calculations dealing with the QD's con- 
tinuum states 2 *^^. This is probably because the non 
trivial shape of the dots renders the continuum states 
hardly analytically solvable. Spherical dot virtual bound 
states were discussed by Buczko and Bassan: 2 while Lc- 
long et a£ related the photoconductivity properties of 
lens-shaped QDs to the existence of virtual bound states. 



In the present paper we examine the structure of 
the electronic continuum states of the InAs/GaAs QDs. 
We shall show that the flatness of the QDs reorga- 
nizes very effectively the continuum states. In fact, it 
is the dominant feature that influences the far infrared 
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This response is associated with 
the photon absorption from the ground state to the con- 
tinuum states. We shall show in particular that a strong 
magnetic field applied parallel to the growth axis leads to 
the formation of quasi-Landau states, as expected. How- 
ever, this quantization is almost invisible in the far in- 
frared response with light polarization along the growth 
direction. Finally, we shall show that one can design 
InAs/GaAs photodetectors which display very strong de- 
pendence of their photoabsorption upon the in-plane po- 
larization of the light. 



II. MODEL 

The InAs/GaAs self-organized dots are known to be 
flat objects (h/R ~ 2-3 nm/10 nm where h is the dot 
height and R the base radius). It is also known that 
they are roughly circular. The energy distance between 
the S'-like symmetry envelope ground state in the con- 
duction band and the center of gravity of the P x and P y 
first excited states is about 50 meV for h about 2 nm 
and R about 10 nmiSiifliil. The splitting between P x 
and P y (about 5 meV) results from potential energy 
terms which do not display cylindrical symmetry, for in- 
stance ellipticity of the QD basis— or/and piezo-electric 
fieldsiii^. The electronic structure of QDs by means of 
the multi-band envelope function method has been cal- 
culated by Stier et ali&ii while accurate atomistic ap- 
proaches (pseudo-potential, tight-binding) were under- 
taken by Williamson and Zunger— (see also Zunger— for 
a review) and by Lee et aUt These numerical meth- 
ods (of either type) mostly aimed at calculating the 
bound states of the dots. Here, we shall use the sim- 
pler (one band) envelope function description to handle 
the more complicated extended states of the dots and 
to predict the QD's far infrared response. We note that 
this one band envelope function method has allowed a 
quantitative description of numerous anti-crossings ob- 
served in the fan charts of the bound-to-bound magneto- 
optical transitions of InAs/GaAs QDs and associated 
with polaronsiSiii. The effective Hamiltonian we inves- 
tigate is: 



H(B) = H (B) + 6V(f) 
P 2 



(1) 

^ + V(p, z) + -oj c L z + ^m*u; 2 p 2 + 5V{r) 

where m* is the carrier effective mass (taken as con- 
stant), V(p, z) the isotropic part of the potential energy, 
lo c = eB/m* the cyclotron frequency and SV the poten- 
tial energy part which does not display the cylindrical 
symmetry. B is the magnitude of the magnetic field ap- 
plied parallel to the z direction. The eigenstates of Ho 
can be chosen as eigenfunctions of L z with the eigenvalue 
to: they will be termed S, P+, P_, D + , £>_, .... A single 
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quantum dot plane faintly absorbs light. Hence, to in- 
crease absorption, one very often uses vertically stacked 
dot planes. When a periodic stacking with a short period 
d (say d = 11 nm) is usedi^, a vertical alignment of the 
QDs is obtained due to their strain distribution 18 . Hence, 
one effectively structures the QD's continuum. This cor- 
responds to looking for Bloch-like solutions of H(B): 

V>fc 2 (s, y,z + d) = e tk * d ■ ^ (x, y, z) (2) 

The eigenenergies in turn become periodic functions 
of k z . The first Brillouin zone will be the segment 
[— n/d, +ir/d]. The beneficial action on the photore- 
sponse of QDs of such vertically periodic stackings with 
the insertion of AlAs barriers in the unit cells has been 
recently investigated^. One possible use of QDs arrays 
is the detection of infrared light. Under such circum- 
stances the dots are modulation-doped and the doping 
concentration is adjusted in such a way that each QD 
contains one electron. The far infrared light is absorbed 
by the QDs because of bound-to-bound transitions (e.g. 
S — > P±) or bound-to-continuum transitions. Two pos- 
sible polarizations may arise: either the electric vector 
£ of the electromagnetic wave is parallel to the z axis 
or £ lays in the layer plane. In the first situation, only 
the continuum states with an S symmetry will be in- 
volved in the absorption processes. In the second situ- 
ation, only the continuum states with a P± symmetry 
will be involved in the response to the electromagnetic 
perturbation (cylindrical symmetry) or, more generally, 
the linear combinations of P + and P_ which diagonal- 
izes H . In the following, we shall present the absorption 
coefficient V versus the photon energy hui for different 
polarizations: 

o * 2 

V{hv) oc V -^^\(S\£- r\v)\ 2 8{e u - e s - hiu) (3) 

where v labels the continuum states. In practice, 
we broadened the delta function by replacing it by a 
Lorentzian of 8 meV full width at half maximum. The 
QD size distribution is known to broaden the absorp- 
tion lines: an average of V over the size distribution of 
the dot (R, h) has to be performed to allow a complete 
comparison to experiments. But the broadening effect 
remains modest for intra-conduction band transitions (a 
few meV) and can easily be taken into accountii, if nec- 
essary. 

A convenient model which takes advantage of the flat 
aspect of the QDs is the separable model^. It will prove 
very useful to interpret the more accurate results ob- 
tained by the numerical diagonalization (see below). In 
the separable model, restricted to the cylindrically sym- 
metrical case for simplicity, the wave function is written 
as: 

An(p,z) = e lm *-g m {p,z) (4) 
g m (p,z) = Ny m (z)f m (p) (5) 



where N is a normalization coefficient, f m (p) & pre- 
scribed function of p which depends on several variational 
parameters (Ai, A2, • . •) and y m {z) an unknown function 
which describes the carrier kinematics along the growth 
axis as resulting from an average potential V m (z): 

V m (z) = N 2 J 2irpdp\f m (p)\ 2 V(p,z) (6) 

V m (z) = V m (z + d) depend on Ai,A2,.... The one 
dimensional Schrodinger equation for the z motion is 
solved numerically. Hence, the eigenvalues of the decou- 
pled problem depends on Ai,A2,..- The lowest bound 
represents the best possible separable solutions. To take 
an example, suppose we discuss the S states: m = 0. We 
try a one parameter Gaussian ansatz for /o: 

fo(p) = exp (-^) (7) 

For a single QD, the lowest eigenvalue (correspond- 
ing to Xmin) will correspond to the ground bound state 
with S symmetry. However, the Schrodinger equation 
for the z motion with V m= ${z, Xmin) admits excited so- 
lutions which can be bound to the QD or belong to the 
continuum spectrum. In the latter case, the decoupling 
procedure amounts to producing a particular set of con- 
tinuum states, those which are bound and nodeless in 
the radial direction but are extended along the growth 
axis. Strictly speaking, these extra solutions should be 
considered with caution because they are not the ground 
state solutions to the problem. However, if the problem 
were truly separable, the variational ansatz would have 
no consequence; the effective potential for the z motion 
would be m independent and all the solutions of the z de- 
pendent problem would be acceptable. Therefore, in the 
case of flat objects such as the QDs, the excited states 
of the effective problem for the z motion may be close to 
actual quantum states of the problem. As we shall see 
in the following they are indeed very useful guidelines. 
Note that in the case of periodic stacking, the excited 
separable solutions will be found to display a significant 
integrated probability of being outside the QD in the unit 
cell of the periodic problem despite the localized nature 
of their lateral motion (as in eq. Q). 

The Gaussian ansatz described above represents the 
15-like states as in atomic physics convention. An or- 
thogonal ensemble of decoupled variational ansatz of 
same (2S, 3S, . . .) or different (1P±, 2P±, . . .) symme- 
try as the ground l^-like one can be easily constructed 
and their eigenenergies are determined following the same 
minimization procedure as for the first one. For instance, 
we chose the following ansatz for the 2S states which have 
one node in the plane: 

/o(p) = (p 2 -^)exp(-^-) (8) 

The po value is related to the parameter X 2 s by imposing 
the orthogonality to the IS states. Hence, the indepen- 
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dent parameter X 2 s was found by minimizing the first 
calculated eigenenergy. 




□ GaAs ^ AlAs □ InAs 



Figure 1: Schematic representation of the supercell including 
the dot and its WL. Period d = 11 nm. 

In the following, we model the QDs as InAs truncated 
cones with basis radius R, height h = 2 nm, a basis 
angle 30° (see fig. QJ. The value of R is taken equal to 
10.2 nm if not mentioned differently. The cones float on a 
two-monolayer (0.5 nm) WL. The dot plane is placed at 
a distance of 1 nm from a 1 nm thick AlAs barrier. The 
overall period is d = 11 nm. This sample was designed 
and investigated experimentally recentl}*^. 

The Bloch eigenstates and eigenvalues of a given m (in 
case SV = 0) are found by projecting the Schrodinger 
equation on a large basis. At zero field, it is a Bessel 
basis for the radial motion^. It corresponds to enclosing 
the QD in a large cylindrical box (100 nm radius) and 
requiring the wave function to vanish at the boundary. 
Plane wave functions with period d are used to describe 
the vertical motion. At high magnetic fields (B > 10 T), 
the Bessel basis is replaced by a set of 2D harmonic os- 
cillator functions (Landau levels). Altogether, we typ- 
ically used 10000 basis functions for the Fourier-Bessel 
basis (21 plane wave functions and 500 Bessel functions) 
and 2800 basis functions for the Fourier-Landau basis (41 
plane wave functions and 70 bidimensional harmonic os- 
cillator functions). The Lanczos algorithm was used to 
extract the 30 lowest lying eigenstates. All calculations 
use a 0.4 eV— conduction band offset between GaAs and 
InAs and 1.08 eV— between AlAs and GaAs. The con- 
duction band effective mass is taken equal to 0.07 mo 
(viz. hu> c = 49.3 meV if B = 30 T) as determined by 
far infrared magneto-absorption experiments^. The en- 
ergy origin is taken at the bottom of the GaAs conduc- 
tion band. We suppose that there is one electron per 
QD coming from doped layers. We also assume that the 
temperature is low enough to take for granted that all 
the electrons are in the ground states of the dots (the 
ground miniband of the stack is dispersionless because of 
the strong localization of the ground state of the individ- 
ual QDs). Consequently, all the absorption processes are 
due to the excitation of these electrons. 



III. RESULTS AND DISCUSSIONS 
A. S states and excitation in z polarization 

Figure |3 shows the calculated levels at k z = and 
k z — ir/d for S states at large field (full lines). We see 
two lowest lying isolated states followed by a series of 
levels organized in fan charts with extrapolations at zero 
field: 16 meV, 144 meV, 293 meV, 701 meV for k z = 
and 21 meV, 101 meV, 371 meV, 560 meV for k z = n/d 
(large dots in fig. |2J). Let us emphasize that these edges 
are roughly those of the GaAs/AlAs/InAs superlattice 
(i.e. with WL but without QDs). With our parameters 
such a superlattice (SL) would display k z — Q (ti / d) edges 
at Si = 19 meV (31 meV), E 2 = 151 meV (110 meV), 
E 3 = 293 meV (374 meV) and E 4 = 699 meV 
(559 meV). The dashed lines in fig. |2 correspond to 
the results of the separable model with a Gaussian ra- 
dial function (eq. Q). The ground state almost coincides 
with the numerical evaluation, as expected. We have 
also checked numerically that the second isolated state is 
indeed a 2S state in the separable model (not shown). 

It is quite remarkable that several extrapolations of the 
fan chart to B = 0, that is to say to the different edges of 
the SL of fig. 2] correspond to the energies of the excited 
separable states with no node in the layer plane [7| Let 
us compare the separable approximation of the actual SL 
in fig. |T]to the GaAs/AlAs/InAs SL (i.e. with WL but 
without QDs). They differ by an extra kinetic energy 
fi 2 1 '(2m* a 2 ) (~ 40 meV) to be paid to localize laterally 
the electron in the dot on the one hand and an extra 
potential energy AV(z) which extends from the WL to 
the top of the dot and increases steadily from -0.4 eV 
to eV— . In a perturbative approach we would have to 
take the average of AV(z) over the probability densities 
associated with wave functions Ck z =a(z) or c ks:=7T u(z) of 
the GaAs/AlAs/InAs SL (see the dashed lines of fig. 
for k z = 0). For the ground state which is nodeless, the 
attractive contribution (ck z =o\ AV\ck z =o) is much greater 
than the lateral kinetic energy. As a result, the first edge 
of the GaAs/AlAs/InAs plunges from +19 meV down 
to -142 meV. Obviously, the first order is insufficient to 
quantitatively account for such a large effect of AV(z). 
In fact we note on fig. [3] that the wave function of the 
separable model is considerably more concentrated on 
the dot along z than the GaAs/AlAs/InAs ground state. 
Going up in energy the second edge (+151 meV) has 
an unperturbed Ck r= o(z) that has a node nearly at the 
center of the period (middle panel in fig. [21 dashed line). 
Hence, the AV attraction is considerably smaller than 
for the ground state, which leads to a much smaller shift 
between the two models. A similar situation holds for 
the excited states. Concomitantly, we see in fig. [3] that 
the wave functions of the separable model become closer 
and closer to those of the GaAs/AlAs/InAs SL. 

The appearance of low lying bound states and high 
energy resonances (dashed lines in fig. |2 see below) 
are an example of the deep attractive perturbation that 
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Figure 3: Probability densities to find the electron at z in 
the unit cell for the three lowest separable solutions with 
a Gaussian radial function and the three lowest states of a 
GaAs/AlAs/InAs SL at B = and k z = 0. From bottom to 
top: the ground state to the second excited state. The verti- 
cal dashed lines delimit the different regions of the structures 
in fig. [U 
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Figure 2: Calculated energy levels of S symmetry states versus 
magnetic field. The dashed lines are the results of the separa- 
ble model with a Gaussian radial function. Upper panel: at 
the center of the first Brillouin zone, lower panel: at an edge 
of the first Brillouin zone. The large dots at B = are the 
extrapolations of the various fan. 



the InAs QD brings to the energy spectrum of the 
GaAs/AlAs/InAs SL. Conversely, the continuum of the 
InAs QD is considerably restructured by the SL effect. 
Because of the small dot height, the resonances (or vir- 
tual bound states) for the z motion occur at very high 
energy. Roughly, one can identify a QD to a quantum 
well with thickness h. Thus, we expect the resonance 
for the z motion to occur at n 2 h 2 /(m*h 2 ) from the InAs 
edge, i. e. in the eV range. Hence, one should expect a 
rather structureless continuum (resonances for the lateral 
motion will be discussed in the next section). Instead, the 
periodicity along z gives rise to an energy spectrum that 
is a periodic function of k z and to wave functions that 
obey the Bloch theorem (eq. EJ). Moreover, the edge of 
the continuum of states is blue-shifted from -15 meV for 
an isolated QD (bottom of the thin WL quantum well) 
to about +16 meV in the periodic structure. The con- 
tinuum states of the QDs and their WLs inserted into a 
GaAs/AlAs SL are therefore unique and cannot be sim- 
ply considered as the result of the perturbation of one of 
the scheme (say the SL) by the other (the QD and its 
WL). 

We show in fig. [21 (full lines) the z variations of the 
probability densities for the three lowest lying separable 
solutions: 



pStUw = \yi:u(z)\ 



(n) 



1,2,3,. 



(9) 



These probability densities are extremely close to 
the actual numerical solutions but have the advan- 
tage of displaying unambiguous nodes while the latter 
is blurred in the exact solutions because of the slight 
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non-separability. The departure from exact separability 
makes that the nodes of the actual problem follow curves 
Tpk z =o{p,z) — 0. When iV'fcJ 2 is integrated over p to 
produce the P^l (z) °f the actual sample, the integra- 
tion over the radial variables changes a node line into a 
minimum along z. The closer this minimum is to zero, 
the more separable the problem is. P^Lo i s nodeless, 

(2) (3) 

PmLo nas one n °de, P^Lq nas two nodes. While the 
ground solution is deeply bound to the dot, the two ex- 
cited solutions leak significantly outside the dot. 

If the separable model were exact, the Landau levels 
which belong to different zone center edges would cross 
at fields such that: 

nhu c +Ei = n'huj c + E v (10) 

where E\ are the eigenenergies for the z motion. 
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have nearly the same in-plane extension. Hence, the re- 
mainder of the continuum states contribute very little to 
the optical absorption. Note that the existence of the 
two peaks in the full calculation and only one in the sep- 
arable model for the first absorption feature at 35 T is 
due to the non-separability of the total wave function, in 
correspondence with the anti-crossing illustrated by large 
circles in fig. [21 Note finally that the absorption spec- 
trum reflects the presence of minibands: their widths 
are equal to E 2 (k z = 0) - E 2 (k z = n/d) « 40 meV 
and E 3 (k z =n/d)- E 3 (k z = 0) « 80 meV for the first 
and second features respectively (the width of the ground 
miniband is negligible) while their double-peak profiles 
reflect the singular density of states at the edges of a SL 
miniband. 
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Figure 4: Absorption coefficients versus photon energy at 
B — and B = 35 T calculated by two models for £ 
parallel to z 



We see in fig. [5] that these crossings are actually re- 
placed by anti-crossings. Their magnitude gives a mea- 
sure of the inaccuracy of the separable ansatz. Since the 
QDs are flat, a great deal of the continuum states can 
be analyzed in terms of separable states. The states pre- 
dicted by the separable model play an important role 
in optical absorption as shown in fig. 0] To obtain 
these results, a summation over 60 k z values of the k z - 
conserving transitions between the initial and final states 
has been performed. The coincidence of the two calcu- 
lated absorption coefficients for the z polarization, using 
either the full 3D calculation or the separable model at 
B = or B = 35 T, evidences that the ground state 
couples preferentially with the continuum states which 



Figure 5: Absorption coefficients versus photon energy from 
B = to B = 35 T for £ parallel to z. 



We show in fig. [SJthe calculated absorption coefficients 
from B = to B — 35 T for light polarized parallel to 
z. The spectrum at B = 10 T is cut off at 0.35 eV due 
to computational limits. We immediately see that the 
spectrum is almost B independent, in spite of a quasi- 
fan chart aspect of the continuum levels (fig. [5J. The 
absorption takes place between the strongly localized ini- 
tial state and the continuum states which 'look' similar 
to the ground state in the layer plane, i. e. the sepa- 
rable excited states (which form a ID continuum with 
frozen in-plane motion and periodic z motion; note that 
the separable model describes excellently the second pro- 
file around 0.45 eV while the marked oscillations in the 
first profile are due to non-separability effects). The tran- 
sitions towards all the other states of the fan charts in 
fig. [21 are very faint because of the mismatched in-plane 
variations between the initial and final states. The in- 
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sensitivity to an external field of the absorption is a di- 
rect consequence of the quasi-separability of the vertical 
and in-plane motions. In this respect the QDs behave 
much the same as narrow one dimensional square wells. 
In these systems, it is known that for £ and B \\ z the 
intersubband absorption is also roughly B independent 
(if it were not for band non parabolicity and/or effec- 
tive mass mismatch there would be an exact decoupling 
between the vertical and in-plane motions and thus an 
exact parallelism between the n th Landau level of all the 
z dependent subbands). 
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B. P states and excitation with in-plane 
polarization 
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Figure 6: Absorption coefficients versus photon energy at 
B — 30 T for £ parallel to z. Usual sample and a centered 
sample (see text). The full line has been rigidly up-shifted for 
clarity. 

Fig. El compares the absorption coefficients for two 
samples: the usual sample has been described above 
while the centered sample has the same material param- 
eters except that the QD has been placed at the center 
of the SL period. By comparing to fig. EJ we see that the 
second absorption peak has disappeared while the first 
one has been reinforced. By placing the QD, which is 
a flat object, at the center of the SL cell, we expect to 
restore a quasi-mirror symmetry with respect to the cen- 
ter of the period, thereby making transitions from the 
ground state (which is even in z) towards even excited 
states almost forbidden in the z polarization. 

To conclude this section, let us emphasize that the ab- 
sorption peaks above are strong and observable in ex- 
periments. At a given magnetic field B and wave vec- 
tor k z , the sum of oscillator strength of the transitions 
from the ground state to the calculated continuum states 
amounts to ~ 0.8. This is close to the saturation value 
(= 1) stated by Thomas-Kuhn-ReichSS sum rule. With 
a structure which has 20 layers and a lateral QD den- 
sity of 5 x 10 10 cm -2 , these peaks were observed either in 
photocurrenlji^ or in absorption* experiments. 



Figure 7: Oscillator strength of transitions from the ground 
S state to the first 30 P states versus energy of P states at 
B — and k z = for several values of base radius R and £ || x. 
The ordinate scale in the case R = 70 A is 5 times bigger 
than the others. The ground energies for each value of R axe: 
E S (R = 70 A) = -105 meV, E S (R = 55 A) = -64 meV, 
E S {R = 50 A) = -45 meV, E S (R = 45 A) = -25 meV 

When the QDs support bound states with P± symme- 
try, the bound-to-bound S — > P± absorption coefficients 
are very strong for £ _L z . These two transitions ex- 
haust almost all the oscillator strength and saturates the 
Thomas-Kuhn-Reich sum rule. Hence only the z polar- 
ization gives rise to a strong bound-to-continuum transi- 
tion (as discussed in section UlI A|l . A possible means to 
overcome the weakness of the absorption with in-plane 
polarization light is to push the P levels in the contin- 
uum where they might form resonant states. Under such 
circumstances the S'-to-continuum absorption will also be 
strong in the in-plane polarizations. This can be achieved 
by decreasing the QD base radius R until the lowest P 
state energy exceeds the onset of the continuum. For 
the ground P level, this happens at about 5.8 nm. Since 
the QD has sharp boundaries, it is likely that the bound 
state will survive in the continuum in the form of a virtual 
bound stat o 3 i 21 . Such a feature would be beneficial to the 
photoconductive properties of the device since, ideally 
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Figure 8: Upper panel: Oscillator strength of transitions from 
the ground S state to the first 30 P states versus energy of 
P states at B = and k z = for 7? = 50 A and £ \\ x. Two 
scales of vertical axis are used below and above 150 meV. The 
apparent scattering of the numerical data above 150 meV is 
due to the contribution of two different optical channels (see 
text). Lower panel: Probability density to find the electron 
at z in the unit cell of some selected states marked (from A 
to E) in the upper panel. 



speaking, a virtual bound state would have a large oscil- 
lator strength distributed over an energy interval broader 
than that of a true bound state and be at the same time 
able to participate in the electrical conduction. To track 
the evolution of the lowest lying continuum P states of 
the QDs, we have calculated for several values of R the 
integrated probability to find the electron in the QD (not 
shown) and the oscillator strength of the S-P transitions 



OSs^p. = 



2m* 



(ep x - e s )\(P x \x\S) | 5 



(11) 



Fig. [7| shows the oscillator strength versus the energy 
of the first 30 final P states for several R values at B = 
and k z = 0. For R = 70 A, one P state is bound. Thus, 
there is a very sharp peak at -14 meV. It is followed by a 
smooth and very weak (OS < 1.5 %) continuum absorp- 
tion in the positive photon energy region. For smaller 
dots, there is no bound P state. Consequently, the oscil- 
lator strength becomes quite strong for final states in the 
continuum spectrum. A peak develops at an increasingly 
larger energy when R decreases. Its width increases with 
decreasing dot size. We note that the sum of the oscilla- 
tor strength inside each broad peak amounts to ~ 0.9, a 
value close to the bound-to-bound transition. This sug- 
gests that the oscillator strength has been redistributed 
over the peak. At a P state energy of about 150 meV, 
independently of R, one notes the existence of two kinds 
of states. A first kind of state appears as the continua- 
tion of the absorption peak while the other kind of state 
are basically optically inactive. We show on fig. [S] the 
evolution of the z dependence for the probability density: 



P r (z) 



2irpdp\ip P!c (p, z)\ 2 



(12) 



with increasing states n. While the states that belong 
to the absorption peak have no node along z, the states 
with a vanishing oscillator strength show one node along 
z. Hence, the apparent cut-off at 150 meV corresponds to 
the onset of a new (but very small) absorption channel, 
this time associated with a change of the quantum num- 
ber of the z motion. We note that the first two eigenen- 
ergies of a SL with GaAs, AlAs, InAs WL but without 
InAs QDs are 19 meV and 151 meV at B = and k z = 
(see discussions for S states in section UlI A|) . that is to 
say the two energies of the absorption onsets shown in 
fig. [7| This shows once again that the flat aspect of the 
QDs leads to a quasi-separability of the carrier motions 
along and perpendicular to the growth axis. 



IV. CONCLUSIONS 

In conclusion, we have calculated the electronic con- 
tinuum states of InAs/GaAs QDs embedded into a 
GaAs/ AlAs SL in the presence of a quantizing magnetic 
field applied parallel to the growth axis. We find that the 
far infrared absorption is essentially independent of the 
vertical magnetic field despite the presence of numerous 
quasi-Landau states. We have proved that the flatness 
of the QDs plays a major part in restructuring the QD's 
continuum. Moreover, we have shown that small dots 
(with a single bound state of S symmetry) may display 
a strong in-plane absorption when a virtual bound state 
with P symmetry is not too far from the continuum edge. 
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Finally, we have shown that most of the optical properties 
of QDs embedded in a SL can be analyzed within a de- 
coupled model. This results from the flatness of the QD 
shape and from the fact that the dipole matrix elements 
are sensitive only to the local (i.e. in the QD region) 
features of the final state wave functions, because of the 
strong localization of the initial QD state. 

Let us finally add a few comments regarding QDs 
with several electrons. It is known that QD can con- 
tain several electrons. When several electrons exist in 
the dots, the Pauli principle may block bound-to-bound 
transitions if the shell of final states is filled. In addi- 
tion, the many electron states in the continuum com- 
prise not only states which can correspond to the de- 
localized states for all electrons but also mixed states 
where some electrons are bound to the dot while the 
others are extended. For instance, in the two electron 
case, the initial state (at low temperature) is a spin sin- 
glet with the two electrons in the ground state. The 
final states can be bound for both electrons or mixed 
with one electron bound and one in the continuum. As 



a result of these mixed states, the (S — S) — (S — P) 
bound-to-bound (discrete) transitions will be accompa- 
nied by (S — S) — (S— continuum) continuum transitions. 
The Coulomb interaction will be more important in the 
bound-to-bound transitions than in the mixed transi- 
tions. Hence, the (S ~ S) — (S — P) discrete lines can 
be very close in energy to the (S — S) — (S— continuum) 
absorption threshold. The description of the absorption 
coefficient of QDs containing more than one electron is 
clearly beyond the scope of this paper. 
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